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ABSTRACT 

The observed samples of supernovae (SN) and double compact objects (DCOs) provide several critical 
constraints on population-synthesis models: the parameters of these models must be carefully chosen 
to reproduce, among other factors, (i) the formation rates of double neutron star (NS-NS) binaries 
and of white dwarf-neutron star (WD-NS) binaries, estimated from binary samples, and (ii) the type 
II and Ib/c supernova rates. Even allowing for extremely conservative accounting of the uncertainties 
in observational and theoretical predictions, we find only a few plausible population synthesis models 
(roughly 9%) are consistent with DCO and SN rates empirically determined from observations. As a 
proof of concept, we describe the information that can be extracted about population synthesis mod- 
els given these observational tests, including surprisingly good agreement with the neutron star kick 
distributions inferred from pulsar proper-motion measurements. In the present study, we find that the 
current observational constraints favor: kicks described by a single Maxwellian with a characteristic 
velocity of about 350km/s (i.e., at maximum likelihood; kick velocities between lOOkm/s and 700km/s 
remain within the 90% confidence interval of unimodal distributions) ; mass- loss fractions during non- 
conservative, but stable, mass transfer episodes of about 90%; and common envelope parameters of 
about 0.15-0.5. Finally, we use the subset of astrophysically consistent models to predict the rates at 
which black hole-neutron star (BH-NS) and NS-NS binaries merge in the Milky Way and the nearby 
Universe, assuming Milky- Way-like galaxies dominate. Inevitably, the resulting probability distribu- 
tions for merger rates depend on our assumed priors for the population model input parameters. In 
this study we adopt relatively conservative priors (flat) for all model parameters covering a rather 
wide range of values. However, as we gain confidence in our knowledge of these inputs, the range of 
merger rates consistent with our knowledge should shift and narrow. 
Subject headings: binaries: close — stars:evolution — stars:neutron - black hole physics 



1. INTRODUCTION 

Many ground-based gravitational wave detectors now 
operating at or near design sensitivity (i.e., LIGO, GEO, 
TAMA, and VIRGO) are designed to detect the late 
stages of double compact object (DCO) inspiral and final 
merger. Based only on early-stage data, these instru- 
ments have already provided conservative upper limits 
to certain DCO merger rates in the Milky Way (see, e.g . 



Abbot et al. (The LIGO Scientific Collaboration! [200l 
Abbott et al. (The LIGO Scientific Collaboration^ 

If LIGO indeed maintains its current level of 
sensitivity (Marx, for the LIGO Scientific Collaboration, 
2006), LIGO data could set rate limits that exclude 
optimistic theoretical predictions, such as a population 
of binary s=s 10M© black holes merging at a present-day 
rate density above ~ lMpc _3 Myr _1 (based on the 
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design LIGO range of 15Mpc for double neutron stars). 
This upper limit has occasionally been reached in 
some theoretical pre dictions for BH-BH m erger rates 
from isolated (e.g. IBelczvnski et al.l I2002L hereafter 
BKB) and interacting (e.g. lO'Shaughnessv et al.ll2007eft 
binary stellar systems. 1 Thus, gravitational-wave based 
upper limits (and, eventually, detections) are expected 
to provide constraints on theoretical models of DCO 
formation. 

Lacking a sample of DCOs containing a black hole, 
the only route to BH-BH merger rates is via popula- 
tion synthesis models. These involve a Monte Carlo 
exploration of the likely life histories of binary stars, 
given statistics governing the initial conditions for bi- 



1 More generally. lO'S haughnessv et al. (2005a c); Brown (1995); 
IPortegies Zwart fc Yungelsonl I|1998I 1: IBelczvnski et al.l l2007ar ); 
and references therein discuss the range of DCO merger rates ex- 
pected to arise from binary evolution in Milky Way-like galaxies. 



2 



naries and a method for following t he behavior of sin- 
gle and binary stars (see, e . g., BKB.lHurlev et alJl2002t 
iPortegies Zwart fc Verbuhtlll996l: [Fryer et al.lll998l and 
references therein). Unfortunately, our understanding of 
the evolution of single and binary stars is incomplete, 
and we parameterize the associated uncertainties with a 
great many parameters (~ 30), many of which can cause 
the predicted DCO merger rates to vary by more than 
one or two orders of magnitude when varied indepen- 
dently through their plausible range. To arrive at more 
definitive answers for DCO merger rates, we must sub- 
stantially reduce our uncertainty in the parameters that 
enter into population synthesis calculations through com- 
parison with observations when possible. 

The simplest and most direct way to constrain the 
parameters of a given population synthesis code is to 
compare several of its many predictions aga i nst ob- 
servations. For example, lO'Shaugh nessv et al.l (|2005d) 
(henceforth PSC1) required their population synthesis 
models be consistent with the empirically estimated for- 
mation rates derived from the three (now four) known 
Galactic NS— NS binaries which are tight enough to 
merge through the emission of gravitational waves within 
10 Gyr. In this paper, we require our population synthe- 
sis models to be consistent (modulo experimental error) 
with six observationally determined rates: (i) the forma- 
tion rate implied by the known Galactic merging NS— NS 
binaries, mentioned above; (ii) the formation rate implied 
by the known Galactic NS— NS binaries which do not 
merge within 10 Gyr (henceforth denoted "wide" NS— NS 
binaries); (iii,iv) the formation rate implied by the sam- 
ple of merging and eccentric WD-NS binaries; and finally 
(v,vi) the type II and type Ib/c SNe rates. Further, we 
use the set of models consistent with these constraints 
to revise our population-synthesis-based expectations for 
various DCO merger rates, assuming no prior informa- 
tion, so all tested population synthesis model parame- 
ters consistent with our constraints are treated equally. 
Though many of the core concepts have been presented 
previously, notably in PSC1, this paper presents more 
constraints, higher-quality fits, and particularly the first 
analy sis of fitting errors (as develo ped in our companion 
paper lO'Shaughnessv et al.ll2007cl henceforth PSF2) to 
arrive at substantially more reliable rate predictions that 
satisfy all available empirical rate constraints. Addition- 
ally, we present both the one and two-dimensional distri- 
butions of those population synthesis parameters which 
reproduce our complete set of constraints. 

We note however, that we focus only on compact ob- 
ject binaries with one or two neutron stars. Merger rate 
predictions for binaries with two black holes must ac- 
count for the mix of both spiral (Milky- Way like) and 
elliptical galaxies in the local Universe, since present- 
day merging double black hole binaries are often sub- 
stantially longer- lived than their lighter counterparts; 
thus, compared to the small fraction of extremely long- 
lived progenitors of merging low-mass binaries, a sig- 
nificant fraction of presently merging BH-BH binaries 
could have been born ~ 10 Gyr ago during the epoch 
of elliptical galaxy formation. As modeling of the lat- 
ter part involves computational elements and method- 
ology that is not part of the present study (i.e., a two- 
component, time-dependent star formation model), we 
intend to undertake BH-BH rate predictions in a forth- 



coming paper lO'Shaughnessv et ail (|2007al ); see also 
Ide Freitas Pacheo et al.1 (|2005| ). However, the evolution- 
ary channels and essential physics ar e very thoroughly 
addressed in IBclczvns ki et al.l (|2007bf ) . 

In Sections [2] and [3] we describe the observational 
constraints we impose on the population synthesis 
models: agreement with empi r ical D CO merger rates 
[NS-NS , based on iKim et all (f2003h (hereafter KKL ) 
and PSC1 ; NS-W D, drawn from iKalogera et all (|2005h : 
Kim et ail (120041) 1 ; and supernova rates [drawn from 
Cappellaro et all (|1999f )]. Considering the significant 
one-sigma systematic error in the pulsar birthrate 
(|Vranesevic et al.1 l2004f ) and the strong correlations of 
that observable to the already included DCO merger rate 
observations, we do not impose the pulsar birthrate as 
an additional constraint. In i|4] we briefly describe the 
version of the Star Track code used and summarize the 
updated results from PSF2 on our population synthesis 
models and the resulting predictions for DCO formation 
rates and SN rates. In Section [5] we select from our fam- 
ily of possible models those predictions which are consis- 
tent with observations of DCOs. We thoroughly review 
the information gained from these constraints and dis- 
cuss their implications for our understanding of binary 
evolution. Also, we employ this sample to generate re- 
fined predictions for the BH-NS and NS— NS merger rates 
through the emission of gravitational waves. 

2. OBSERVATIONS OF NS DCOS 

Seven NS— NS binaries and four WD-NS binaries 
(with relatively massive WDs) have been discovered 
so far in the galactic disk, using pulsar surveys with 
well-understood selection effects. We are very specifi- 
cally not inc luding the recently- discovered binary PSR 
J1906+0746 (|L primer et al.l 1200(f ), because the compan- 
ion cannot be definitively classified as a WD or NS at 
present 2 . We also omit PSR J2127+11C found in the 
globular cluster M15, since its formation is thought to 
be dynamical and not due t o isolated binary evolu tion in 
the Galactic field; (see,e,g. lAnderson et al.l fl990). The 
basic measured and derived properties of these binaries 
are presented in Table [TJ 

2.1. Methodology I: Preferred population model 

In the context of NS— NS , KKL developed a statisti- 
cal method to estimate the likelihood of DCO formation 
rates, given an observed sample of DCOs in which one 
member is a pulsar, designed to account for the small 
number of known systems and their associated uncertain- 
ties. Their analysis assumes that each binary is a unique 
representative of its own class. For each primitive class 
v (i.e., for each pulsar), they find a posterior probability 
distribution function V v for its formation rate 1Z U : 

V v {1l)=Allle- A " n . (1) 

The parameter A v depends on some of the properties 
of the pulsars in the observed DCO sample [see KKL 
Eq. (17)]: 

A = m c /(f b N PSR ) (2) 

2 Our predictions for the galactic NS-NS merger rate would in- 
crease by a factor of roughly two if indeed PSR J1906+0746 were 
a double neutron star. 
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Table 1. Properties of compact binaries 



PSRs" P s b P s c P b d Mf, SR u( c» /* r* r J d r* rg rjj N™ SR A" Refs° 
H (10- 18 ) (hr) (Mq) (M q ) (Gyr) (Gyr) (Gyr) (Gyr) (Myr) 


(1) NS-NS(m) 




B1913+16 59.03 8.63 7.752 1.4408 1.3873 0.617 5.72 
B1534+12 37.90 2.43 10.098 1.3332 1.3452 0.274 6.04 
J0737-3039A 22.70 1.74 2.454 1.337 1.25 0.088 6* 
J1756-2251 28.46 1.02 7.67 1.4 1.18 0.181 6* 


0.11 0.065 0.3 4.34 620 0.103 1,2 
0.25 0.19 2.7 9.55 440 1.011 3,4,5 
0.207 0.145 0.085 14.1 1620 0.023 6,7 
0.44 0.38 1.7 16 - - 8 


(2) NS-NS(vw) 




J1811-1736 104.182 0.916 451.20 1.62 1.11 0.828 6* 
J1518+4904 40.935 0.02 207.216 1.56 1.05 0.25 6* 
J1829+2456 41.0098 ~ 0.05 28.0 1.14 1.36 0.139 6* 


1.8 1.74 n/a 7.8 606 2.61 9,10 
32 32 n/a 54.2 282 37.8 11 
13.0 12.9 n/a 43.3 272 32.7 12 


(3) WD-NS(m) 




J0751+1807 3.479 0.00779 6.312 2.1 0.191 ^0 6* 
J1757-5322 8.870 0.0263 10.879 1.35 0.67 « 6* 
J1141-6545 393.9 4290 4.744 1.30 0.986 0.172 6* 


7.08 6.9 4.3 400 2850 0.66 13,14 
5.06 4.93 8.0 145 1240 1.74 15 
0.0015 - 0.56 0.104 370 0.047 16,17 


(4) WD-NS(e) 




J1141-6545 393.9 4290 4.744 1.30 0.986 0.172 6* 
B2303+46 1066.37 569 295.2 1.34 1.3 0.658 6* 


0.0015 - 0.56 0.104 370 0.047 15,16 
0.0297 0.014 n/a 0.138 240 0.106 17,18 



"Name of double-compact binary observed via a component pulsar. 
'Spin period. 
c Spin-down rate. 
d Orbital period. 

^Estimated mass of visible pulsar. 
^Estimated mass of the pulsar companion. 
^Eccentricity. 

^Beaming correction factor. An asterisk indicates the beaming factor has not been experimentally determined; for these systems, we use 
the canonical value /(, = 6 (shown). 
'Characteristic age of the pulsar. 

^Spin-down age of the pulsar. We calculate T sd only for recycled pulsars. 
^Merging time of a binary system due to the emission of gravitational waves. 
'Death time of the pulsar. 

m Most probable value of the total number of pulsars in a model galaxy estimated for the reference model (model 6 in KKL). 

"Parameter in rate equat ion [see Eq, j2H . 

" References: ( 11 IHulse fc Tavlorl Hl9W): ( 2) IWex et al.l iHfflOOTl ; (31 IWolszczanl fl99lT ) : (41 IStairs et all (f2003): ( 51 IStairs et alj Pool): 

(61 ILvne et al. ( 20041 ~(7llWolsz czan ( 1991 ]); (8llFaulkner et al.) I|2005I1; (9) ILvne etaJ] (I2000l ) ; (10) IKramer et al.| (120031 ); fil l INice et al. 

1 1996 1 ) ; (12) ICharnpiori ot alj ^0041) 7(131 Lund gren et alj (11995(1; (141 INice et alj (120051) ; (15) Edwards & Bailcs (2001); (16) Kas pi et al 

(2000): (17) lBailes et alj 1 12003 ); (18l lStokes et al.l (I1985T) ; (19) Ivan Kerkwiik fc Ku lkarni (199g); 



where fr is the fraction of all solid angle the pulsar 
beam subtends; riif c is the total detectable lifetime of 
the object; and JVpsr is the total estimated number of 
systems similar to each of the observed ones (i.e., Np§ K 
is effectively a volume-weighted probability that a pul- 
sar with the same orbit and an optimally oriented beam 
would be seen with any of the radio pulsar surveys; this 
factor incorporates all our knowledge of pulsar survey se- 
lection effects as well as the pulsar space and luminosity 
distributions). In this section, we present results only 
for our preferred pulsar luminosity distribution model, 
corresponding a power law luminosity distribution with 
negative slope, index p = 2, and minimum "luminosity 3 " 
Lmin = 0.3 mJy kpc 2 [model 6 of KKL; as discussed 
therein, this model better accounts for more recent ob- 
servations of faint pulsars]. In the following section, we 
describe how our predictions change when systematic un- 
certainties in p and L m i n are incorporated into V '. 

We obtain -/Vpsr via detailed simulations of pulsar 
survey sel ection effects, as des cribed in KKL (see also 
KKL and Kalogcr a et al.ll2005f ). We obtain the beam- 
ing correction factor via direct pulsar measurements 
(jKalogera et al.ll2001l ) for those systems for which data is 
available. In the cases where /& is not known, we choose a 

3 More correctly, this quantitiy represents pseudo-luminosity us- 
ing a 400 Mhz observing frequency. 



canonical value fb = 6 consistent with existing measure- 
ments. Finally, we estimate the detectable lifetime riif e 
via pulsar and binary properties. The relevant lifetime 
for binaries detected via a pulsar component is 

= min(r sd , lOGyr) + min(r d , T mrg ) (3) 

(where T s d is t he p ulsar spindown age [see 
I Arzoumanian et "ail (|1999[ )]; T mrg is the time remaining 
until the binary merges through the emission of gravita- 
tional waves [see iPetersI (|1964l) and iPeters fc Mathewi 
(119631)1; and Td is the puls ar "death time" [Eq. (9) from 
IChen &; Rudermanl (|1993l )]l. This calculation assumes 
the binary will be detectable via pulsar emission from its 
birth (estimated by the spindown contribution) up until 
its eventual merger via the emission of gravitational 
waves (estimated via T mrg ) or until the pulsar ceases 
to emit (estimated by Td). [Since r s d estimates are 
somewhat uncertain, we require that they do not exceed 
the current age of the Galactic disk (10 Gyr).] Table [1] 
presents appropriate values for /;,, iVpsR, Thfc an d A for 
each individual pulsar-containing DCO. 

Using these individual rate estimates, we then generate 
a distribution for each net formation rate for any class K 
by assuming our individual representatives are complete 
as well as independent - in other words, we assume any 
any member of class K (for K =NS-NS, any merging 
double NS) is similar to one of the observed members 
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Table 2. Constraint intervals 



Type 


Low 


High 




Pok+ 




1 °gio( i? y r ) 


1 °gio(- R y r ) 




NS-NS(m) 


-4.50 


-3.65 


0.26 


0.40 


NS-NS (vw) 


-6.62 


-5.71 


0.36 


0.60 


WD-NS(e) 


-4.94 


-3.98 


0.45 


0.65 


WD-NS(m) 


-4.70 


-3.91 


0.35 


0.60 


SN Ib/c 


-3.34 


-1.97 


1 


1 


SN II 


-2.60 


-1.55 


1 


1 



"Fraction of population synthesis model parameters in our canon- 
ical model space which produce rates in the constraint-satisfying 
interval. See § [4] 

"Fraction of population synthesis model parameters which pro- 
duce rates within the constraint-satisfying interval, modulo an es- 
timate of the rms systematic error associated with our fitting pro- 
cess. See §H and PSF2. 

of that class. [Lacking definitive information about or- 
bital and pulsar parameter distributions, we prefer not to 
perform unwarranted extrapolations to account for bina- 
ries that have not been seen.] For example, to estimate 
the formation rate of merging NS— NS binaries in Milky 
Way- equivalent galaxies (MWEG) , we take three merg- 
ing NS— NS binaries and combine their formation rate 
via 

T'ns-ns (Ktot) = J dn x dTl 3 dTl 2 ^Tl tot -n x -Tl 2 - TZ 3 ) 

xPi(Ki)P2(K 2 )P 3 (K 3 ) (4) 

[This expression is discussed in §5.2 of KKL and pre- 
sented exp l icitly for the three-binary case in Eq. (A8) of 
iKim et all l|200l .] 

Finally, for each class K and thus each class rate distri- 
bution Vk we define symmetric 95% confidence intervals: 
the upper and lower rate limits ~R W ,± satisfy 

/ dKV w (K)= dVJP w {K) = 0.025 . (5) 

Jo Jn Wi+ 

This confidence-interval convention is different from the 
customary choice presented in KKL. 

In all plots that follow, rate probability distributions 
are represented using a logarithmic scale for TZ; thus, 
instead of plotting P, all plots instead show 

p(logft) = V{1Z)11\q.10 . 

2.2. WD-NS binaries 

Four WD-NS binaries with relatively massive WDs 
have been discovered in the galactic disk: PSRs 
J0751+1807, J1757-5322, Jl 141-6545, and B2303+46. 
The intrinsic and deduced properties of these four 
WD— NS binaries are summarized in Table Q] While 
all four binaries may be applied to a net WD-NS for- 
mation rate estimate, the sample manifestly contains 
relics of distinctly different evolutionary channels; for 
example, while J0751+1807 and J1757-5332 have evi- 
dently been strongly circularized and spun up by a re- 
cent mass transfer e pisode, Jl 141 -6545 cannot have been 
(jKaspi et alj 120001 : Bailes et all I2003D . Ideally, popu- 
lation synthesis must produce distributions of WD-NS 
binaries consistent with both the observed orbital pa- 
rameters and spins. The present sparse sample, how- 
ever, does not allow a reliable nonparametric estimate 



of the distribution of WD-NS binary parameters. In- 
stead, as first step towards applying constraints based 
on binary parameter distributions ; we subdivide these 
four binaries into two overlapping classes: merging bi- 
naries, denoted WD-NS(m), which will merge through 
the emission of gravitational waves within 10 Gyr; and 
eccentric binaries, denoted WD-NS (e), which have sig- 
nificant (e > 0.1) eccentricity at present. The rate esti- 
mate derived for both classes is domina ted by J1141-6545 
(|Kim et al . 2004: Kalog era, et al.1 2005) . Ignoring system- 
atic uncertainties (e.g., assuming the beaming correction 
factor ft, is indeed its canonical value and not account- 
ing the pulsar luminosity function uncertainties), Tabled 
shows the 95% confidence interval for each class' forma- 
tion rate, shown per year for a Milky Way-like galaxy. 

2.3. NS-NS binaries 

Seven NS— NS binaries have been discovered so far 
in the Galactic disk. Four of the known systems will 
have merged within 10 Gyr (i.e., "merging" binaries: 
PSRs J0737-3039A, B1913+16, B1534+12, and J1756- 
2251) and three are wide with much longer merger 
times (PSRs J1811-1736, J1518+4904, and J1829+2456). 
The intrinsic and deduced properties of these seven 
NS-NS binaries are summarized in TableQ] PSR J1756- 
2251 was discovered recently with acceleration searches 
(jFaulkner et al.ll2005fh Since its pulsar and orbital pa- 
rameters are sufficiently similar to PSR 131913+16 and 
since the selection bias due to acceleration was already 
included in KKL, we do not expect this new system to 
significantly affect the NS-N S inspiral rate esti mates. As 
discussed in more detail by IKim et ail ((2006) , the rate 
increase is estimated to be smaller than 4%. For this 
reason, we omit it when estimating NS-NS merger rates. 

The observed NS-NS population naturally subdivides 
into two distinct classes, depending on whether they 
merge due to the emission of gravitational waves within 
10 Gyr: merging NS-NS binaries, denoted NSNS(m), 
and wide NS-NS binaries, denoted NSNS(vw). Table [2 
shows the 95% confidence interval for each class' forma- 
tion rate. 

2.4. Recycling, selection effects, and the lack of wide 
NS-NS binaries 

As Table [2] indicates, the confidence intervals for wide 
and merging NS-NS binaries are almost an order of mag- 
nitude from overlapping. On the contrary, population 
synthesis simulations produce merging and wide bina- 
ries at a roughly equal (and always highly correlated) 
rates. Since our rate estimation technique automatically 
compensates for the most obvious selection effects (e.g., 
orientation and detectable lifetime), two unbiased sam- 
ples of wide and merging NS-NS binaries should arrive 
at nearly the same prediction for the NS-NS formation 
rate. 

PSC1 explained this significant discrepancy by a selec- 
tion effect: evolutionary tracks leading to wide NS-NS 
binaries should be less likely to recycle the first-born 
pulsar. The observed sample of wide NS-NS binaries 
[NSNS(vw)] represents a much smaller subset of recy- 
cled pulsars. As summarized in PSF2, we assume any 
wide NS-NS binary which in its past underwent any con- 
ventional, i.e., dynamically stable, mass transfer will re- 
cycle one of its neutron stars to a long-lived pulsar. This 
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condition is likely to over-estimate the true NSNS(vw) 
formation rate. We note that neutron-star recycling can 
also possibly happen during dynamically unstable mass- 
transfer phases (common envelopes), but the vast major- 
ity of wide NS-NS do not evolve through such a phase. 

2.5. Pulsar Population Model Uncertainties 

As noted in KKL and subsequent papers, our recon- 
struction of the pulsar population (i.e., Npsr) relies 
upon our understanding of pulsar survey selection ef- 
fects and thus on the underlying pulsar luminosity dis- 
tribution. This distrib ution can be well-constr ained ex- 
perimentally (see, e.g. lCordes fc Chernofj|l997t ). though 
these constraints do not yet incorporate re cent faint pul- 
sar d iscoveries such as PSR J1124-5916 ()Camilo et al.l 
2002). Nonetheless, different observationally-consistent 
distributions imply significantly different distributions, 
with maximum-likelihood rates differing by factors of or- 
der 10 (KKL). Since the constraint intervals discussed 
above assume the preferred pulsar luminosity distribu- 
tion model - a power-law pulsar luminosity distribution 
with negative slope p = 2 and minimum cutoff luminos- 
ity L m i n — 0.3 mjy kpc 2 - they do not incorporate any 
uncertainty in the pulsar luminosity function. 

The infrastructure needed to incorporate uncertainties 
in the pulsar luminosit y function has bee n presented and 
applied, for example in lKim et al.l (|2006f ). However, out- 
of-date constraints on the pulsar luminosity function al- 
low implausibly high minimum pulsar luminosities L m .;„. 
A high minimum pulsar luminosity implies fewer merging 
pulsars have been missed by surveys. Thus, these out- 
of-date constraints on pulsar luminosity functions permit 
models consistent with substantially lower merger rates 
than now seem likely, given the discovery of faint pulsars. 
In other words, i f we use the infrastructure presented in 
iKim et all (|2006f) to marginalize over Lmin and p gener- 
ate a net distribution function for the merger rate, then 
the 95% confidence intervals associated with that net dis- 
tribution would have a spuriously small lower bound, en- 
tirely because the pulsar population model permits large 
Lmin- Therefore, in the present study we present re- 
sults based only on our preferred luminosity function and 
do not inc lude the out-of-date lumin osity function con- 
straints bv lCordes fc Chernoffl (119971) and the r elated net 
rate distribution provided by IKim et al.l (|2006| ). 

3. OBSERVATIONS OF SUPERNOVAE 

Type Ib/c and II supernovae occur extremely rarely 
near the Milky Way. While historical data contains sev- 
eral observations of and even surveys for supernovae, 
the selection effects in these long-duration heteroge- 
neous data sets have made their interpretation difficult 
(Cappellaro 2005, private communication). In this pa- 
per, we estimat e supernova rates and u ncertainties us- 
ing Table 4 of ICappellaro et al.l (|1999h . To translate 
this table, which expresses results and their (one-sigma, 
systematics-dominated) uncertainties in terms of a num- 
ber of supernovae per century per 10 10 blue solar lu- 
minosities (Lq^b), into the equivalent rates per Milky 
Way equivalent galaxy, we assume a Milky Way blue lu- 
minosity of 1.7 x 1O 1o L 0iS = 0.9 x 10 10 L© relative to 
the solar blue (Lf^_p) and total (L^ ) lumi nosities (see 
e.g., 



pernovae rates that appear in their Table 4, expressed 
in "SNU" (1 SNU = l/10 10 L QjB x 10 2 yr) into a super- 
novae rate expressed in events per year per Milky Way- 
equivalent galaxy (MWEG) using a proportionality con- 
stant W = 2 x 10- 2 MWEG/SNU. For simplicity, we 
ignore O(30%) uncertainty i n the blue luminosity of the 
Milky Way and thus W (see iKalogera et al1l200ll for an 
estimate of this uncertainty) . 

Unfortunately, without substantially more detailed in- 
formation on and statistics of the limiting (systematic) 
errors, we cannot justify any relation between the lin- 
ear, one-sigma confidence intervals provided in their Ta- 
ble 4 and the logarithmic, two-sigma confidence intervals 
needed in our analysis. Lacking any way to rigorously 
justify a useful choice 4 for a bounded two-sided confi- 
dence interval in logi?, we decided to employ the loga- 
rithmic confidence interval 



log(rVF) ± 2 max 



log 



r + u 



r 



log- 



(6) 



juiai ume \^(7) hi aim uuuai yijp) ) lummosiues ^see, 
IKalogera et aHl(2001l ). [Phinnevl l| 19911 ), ?, and ref- 



erences therein). More specifically, we translate the su- 



where r and u are the average value of and one-sigma 
uncertainty in the supernovae rate (in SNU); the corre- 
sponding explicit confidence intervals (in logi? per year 
per MWEG) for Type II and Ib/c supernovae are shown 
both in Table and Figure [TJ Intuitively, this logarith- 
mic interval arises from (i) converting the mean value and 
limiting bounds of a one-sigma confidence interval in r to 
a "one-sigma" confidence interval in logi?; (ii) calculat- 
ing the maximum distance from the central value to the 
corresponding two "one-sigma" limits in log i?; and then 
(iii) doubling that distance to generate a "two-sigma" 
interval. [A marginally more optimistic case using the 
minimum distance from center to each limit also yields 
nearly the same intervals.] 

Though fairly accurate stu dies exist of the high- 
redshift supernova rate (e.g., ICappellaro et al.l 120051) . 
they have little relevance to the present-day Milky Way. 
Several surveys have also attempted to determine the su- 
pernova rate in the Milky Way by a variety of indirect 
methods, such as statistics of supernova remnants (highly 
unreliable due to challenging selec tion effects; see, e.g., 
Ivan den Bergh fc Tammannl Il99lh and direct observa- 
tion of radioisotope-produ ced backgrounds (e.g., decay 
from 26 Al, as described in Die hl et aLll2006f ). Taken in- 
dependently, these methods have greater uncertain ties 
than the historical studies of ICappellaro et ah! ([T999). 

4. POPULATION SYNTHESIS PREDICTIONS 

4.1. StarTrack population synthesis code 

We estimate formation and merger rates for sev- 
eral classes of double compact objects using the Star- 
Track code first developed by BKB and recently sig- 
ni ficantly updated and te sted as de scribed in detail 
in iBelczvnski et all (|2006bf ); see §2 of IBelczvnsk i et all 
(|2007b| ) for a succinct description of the changes between 
versions. This updated code predicts somewhat different 

4 In fact, the Chebyshev inequality P(\x — fi\ > t) < <r 2 /t 2 , which 
applies to any random variable with mean fi and standard deviation 
a, provides a universal tool with which we can find a better-than- 
95% "overconfidence interval" on r and thus logi?. However, not 
only does that interval in r includes r = - so arbitrarily small 
supernovae rates would be included - but it also corresponds to 
an excessively conservative 4.5cr confidence interval for Gaussian- 
distributed errors. 
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double compact object properties than the version used 
in BKB. 

Additionally, we have adopted a maximum neutron 
star mass of 2.5M Q (raised fro m 2M& used in our 
two pr evious studies PSC1 and lO'Shaughnessv et al.l 
(|2005af ). henceforth PSF1), motivated by recent obser- 
vati ons sugges t ing pu lsars with masses near 2M P) (see, 
e.g. lNice et all (|2005t ) and ISplaver et all (|2002D 1. Such 
a higher maximum neutron star mass effectively con- 
verts many merging binaries we would otherwise in- 
terpret as BH-NS or even BH-BH binaries into merg- 
ing NS binaries. Consequently the BH-NS and partic- 
ularly BH-BH rates are decreased from the values seen 
in PSC1. [As a concrete example, with a low maximum 
NS mass, many BH-NS binaries form through accretion- 
indu ced collapse of a NS to a B H, as seen in Figure 
1 of lO'Shaughnessv et alj (|2005bD . However, a signifi- 
cant fraction of these systems would be interpreted as 
NS-NS bi naries with the adopted higher maximum NS 
mass; see [Bclczyns ki et al.1 (|2007bl ) for a detailed discus- 
sion to merger rates that this revised minimum neutron 
star mass implies.] 

Like any population synthesis code, it evolves ran- 
domly chosen binaries from their birth to the present, 
tracking the stellar and binary parameters. For any class 
of events that is identifiable within the code, such as su- 
pernovae or DCO mergers, we estimate event rates by 
taking the average event rate within the simulation (i.e., 
by dividing the total number of events seen within some 
simulation by the duration of that simulation) and renor- 
malizing by a scale factor that depends on properties 
of the simulation (i.e., the number of binaries simulated 
and the binary birth mass distributions assumed) and 
the Milky Way as a whole (i.e., the present-day star for- 
mation rate). Specifically, our simulations are normal- 
ized to be consistent with an assumed Milky Way star 
formation rate 5 of 3.5M0yr _1 . The necessary formulae 
are presented and discussed extensively in a companion 
paper PSF2. 

Finally, like any population synthesis code, the rela- 
tive probabilities of different evolutionary endpoints are 
heavily influenced by the relative likelihood of different 
initial conditions, such as the distribution of initial semi- 
major axes. For example, observational data is consis- 
tent with initial semimajor a xes distribu ted uniformly 
in log q; see e.g., Figure 2 oflAbtJ |l983) and Figure 7 
of iDuauennov fc Mavorl (|1991l ). 6 Considering the limits 
of the period distribution and the masses involved, we 
adopt a semimajor axis distribution flat in logo from 
contact to 10 5 i? Q (see BKB). 

When constructing our archived population synthesis 

5 Our approach gives only the average event rate. The present- 
day merger rate agrees with this quantity when most mergers oc- 
cur relatively promptly (e.g., < 100 Myr) after their birth. Some 
DCOs - notably double BH binaries - have substantial delays be- 
tween birth and merger, introducing a strong time dependence to 
the merger rate. The technique described above will significantly 
underestimate these rates. This point will be addressed in consid- 
erably more detail, both for the Milky Way a nd for a heterogeneous 
galaxy population in forthcoming papers by O'Shaughncssy et al. 
(2007111). 

6 The semimajor axis distributions in these two papers are rela- 
tively similar, though they do predict a slightly different proportion 
of binaries in any small bin dn/da. However, for spiral galaxies the 
instantaneous merger rate depends on the fraction of binaries inside 
a broad period interval, not the derivative at any point. 



results, we did not choose to record detailed information 
about the nature and amount of any mass transfer onto 
the first-born NS. We therefore cannot reconstruct pre- 
cise population synthesis predictions for the NSNS(vw) 
formation rate. However, we do record whether some 
mass transfer occurs, and the nature of the mass transfer 
mechanism. The conventional mass transfer mechanism 
- dynamically stable Roche-lobe overflow - inevitably 
produces a disk around the compact object and can po- 
tentially spin that object up. Other mechanisms, such 
as (possibly hypercritical) common-envelope evolution, 
presumably involve a substantially more spherical accre- 
tion flow; the specific angular momentum accreted may 
be substantially lower, possibly not enough to recycle the 
neutron star. Thus for the purposes of identifying a class 
of potentially recycled ("visible") wide NS-NS binaries, 
we assume any system which underwent mass transfer 
(dynamically stable in the case of wide systems) pro- 
duced a recycled NS primary. 

4.2. Fitting results 

As in previous studies such as PSF1, here we have cho- 
sen to vary seven (7) model parameters in the synthesis 
calculations. These choices are strongly guided by our 
past experience with double-compact-object population 
synthesis and represent the model parameters for which 
strong dependence has been confirmed. These seven pa- 
rameters enter into every aspect of our population syn- 
thesis model. One parameter, a power law index r in 
our parameterization of the companion mass distribu- 
tion, influences the initial binary parameter distributions 
(through the companion masses). Another parameter w, 
the massive stellar wind strength, controls how rapidly 
the massive progenitors of compact objects lose mass; 
this parameter strongly affects the final compact-object 
mass distribution. Three parameters U1.2 and s are used 
to parameterize the supernovae kick distribution as a su- 
perposition of two independent maxwellians. These kicks 
provide critical opportunities to push distant stars into 
tight orbits and also to disrupt potential double neutron 
star progenitors. Finally, two parameters a A and f a gov- 
ern energy and mass transfer during certain types of bi- 
nary interactions; see Section 2.2.4 of BKB for details. 

Other parameters, such as the fraction of stellar sys- 
tems which are binary (here, we assume a binary fraction 
of 100%) and the distribution of initial binary pa ram- 
eters, ar e comparatively well- d etermi ned (see e.g. lAbtj 
(1983), IDuauennov fc Mavorl (|1991l ) and references 
therein); for these parameters we have adopted preferred 
values. 7 Even for the less-well-constrained parameters, 
some inferences have been drawn from previous studies, 
either more or less directly (e.g., via observations of pul- 
sar proper motions, which presum ably relate closely t o 
supernovae kick stre n gths; see , e.g., iHobbs et al. | (2 005). 
lArzoumanian et ail (|2002D , iFaucher-Giguere fc Kaspil 
(2006) and references therein) or via comparison of 
some subset of binary popu lation synthesi s result s 
with observations ( e. g., §8 o f iBelczvnski et all (12006b). 
Ivan der Sluvs et all (|2006l) . INelemans fc ToutI (|2005l) . 

7 Particularly for the application at hand - the gravitational- 
wave-dominated delay between binary birth and merger - the de- 
tails of the semimajor axis distribution matter lit t le. For a similar 
but more extreme case, see lO'Shaughnessv et al.l (l2007elV 
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iWillems fc KolTJ (pOO^ . IPodsiadlowski et all (pOOl) and 
references therein). Despite observational suggestions 
that point towards preferred values for these seven pa- 
rameters - and despite the good agreement with short 
GRB and other observations obta i ned wh en using these 
preferred values (jBelczvnski et alJ (|2006aT) ) - in this pa- 
per we will conservatively examine the implications of 
a plausible range of each of these parameters. More 
specifically, despite the Milky Way -specific studies of 
lO'Shaughnessv et all (|2005d l2007dl ) (which apply only 
to spirals, not the elliptical galaxies included in this pa- 
per), to allow for a broad range of possible models, we 
consider the specific parame ter ranges quoted in §2 of 
lO'Shaughnessv et al.l (|2005al ) for all parameters except 
kicks. For supernova kick parameters, we allow the dis- 
persion v\ t 2 of either component of a bimodal Maxwellian 
to run from to 1000 km/s. 

Since population synthesis calculations involve con- 
siderable computational expense, in practice we esti- 
mate the merger rate we expect for a given combi- 
nation of population synthesis model parameters via 
seven-dimensional fits to an archive of roughly 3000 de- 
tailed simulations, as presented and analyzed in detail 
in PSF2. However, these fits introduce systematic errors, 
which have the potential to significantly change the pre- 
dicted set of constraint-satisfying models. For this rea- 
son, PSF2 also estimates the rms error associated with 
each seven-dimensional fit. Finally, PSF2 demonstrates 
that, by broadening the constraint-satisfying interval 
by the fit's rms error, the predicted set of constraint- 
satisfying models includes most models which actually 
satisfy the constraints. For this reason, the dark shaded 
regions in Figures [TJ [2j and [3] account for both observa- 
tional uncertainty in the Milky Way merger rate and for 
systematic errors in the fits against which these observa- 
tions will be compared. 

4.3. Prior distributions 

Lacking knowledge about which population synthe- 
sis model is correct, we assume all population synthe- 
sis model parameters in our range are equally likely. A 
monte carlo over the seven-dimensional parameter space 
allows us then to estimate the relative likelihood, all 
things being equal, of various DCO merger rates (shown 
in Figures 1 21: j L and [6]) and supernova rates (Figure [lj . 
Also shown in these figures are the observational con- 
straints described in Sections and [3] (shown in shaded 
gray). Table [2] provides both the explicit confidence in- 
tervals of the constraints and the fraction of a priori 
models which satisfy that constraint. Finally, the dotted 
lines in Figure [6] show the a priori population synthesis 
predictions for BH-NS and NS-NS merger rates in the 
Milky Way. 

5. APPLYING AND EMPLOYING CONSTRAINTS 

Physically consistent models must reproduce all pre- 
dictions. Supernovae rates, being ill-determined due to 
the large observational systematic errors mentioned in 
Sec. [3j are easily reproduced by almost all models at the 
2cr level; see Fig. [TJ However, population synthesis mod- 
els do not always reproduce observed formation rates for 
DCOs, as shown in Figs. [2] and [3j These four 95% con- 
fidence intervals implicitly define a (0.95) 4 w 81% confi- 
dence interval in the seven-dimensional space, consisting 
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Fig. 1. — Probability distributions for SN Ib/c (top) and SN 
II (bottom), as predicted by population syntheses with StarTrack. 
The light shaded region shows the (as 2a) interv al consistent with 
observ ational constraints of supernovae, from Cappcllaro et al. 
(1999). A dark shaded region (barely visible in this plot) indicates 
these constraints, augmented by an estimate of the systematic er- 
rors in our fits; see PSF2. These constraints provide no information 
about population synthesis. 



of 9% of the original parameter volume when systematic 
errors in our fitting procedure are included; see Tabled 8 
In other words, we are quite confident (80% chance) that 
the physically-appropriate parameters entering into Star- 
Track can be confined within a small seven-dimensional 
volume, in principle significantly reducing our model un- 
certainty. 

In Figure [4] we show one-dimensional projections onto 
each coordinate axis of the constraint-satisfying volume 
- in other words, the distribution of values each indi- 
vidual population synthesis parameter can take. For 
the purposes of this and subsequent plots, we use the 
following seven dimensionlcss parameters Xk that run 
from to 1: x\ = r/3; x 2 = w; x 3 — ui/(1000km/s), 
X4 = U2/(1000km/s); x$ — s; xq — aA, and 27 = f a . 
As seen in the top panel of this figure, the single most 
likely distribution of supernovae kicks bears a surpris- 
ing rese mblance to the observed pu ls ar kick distribution 
Csee.e.g. lArzoumanian et all 119991 : iHobbs et all 120051 : 
iFaucher-Giguere fc Kaspill2006t and references therein), 
even though no information about pulsar motions have 
been included among our constraints and priors. How- 
ever, as is clear from Figure [T] and more clearly by the 

8 Our confidence in this interval is significantly higher than 81%, 
because we have very conservatively incorporated estimates for sys- 
tematic errors in our fits; see PSF2 for details. 
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Fig. 2. — Probability distributions for merging (top) and eccen- 
tric (bottom) WD-NS formation rates, as predicted by population 
synthesis. The light shaded region shows the 95% confidence in- 
terval consistent with observations of binary pulsars, as described 
in Sec. [2] The dark shaded region extends this constraint interval 
by an estimate of the systematic error in each fit; see PSF2. 



top right panel of Figure [5l a wide range of kick ve- 
locity distributions remain consistent with observations, 
including those with very low (~ lOOkm/s) and very 
high (~ 600km/s) characteristic velocities. The popula- 
tion synthesis mass transfer parameter f a is particularly 
well-constrained, with a strong maximum near x = 0, 
implying that mass-loss fractions of about 90% or higher 
are favored in non-conservative, but stable, mass transfer 
episodes. Also, common envelope efficiencies aX of about 
0.15-0.5 appear to be favored by the constraints. The 
rest of the one-dimensional parameter distributions are 
nearly constant and thus uninformative. Though small, 
the constraint-satisfying volume extends throughout the 
seven-dimensional parameter space. Higher-dimensional 
correlations disguise most of the remaining information, 
as indicated in Figure [5] 

These results should be contrasted with the dis- 
tinctly different post-cons traint results shown in 
lO'Shaughnessv et all (|2005cD (cf. their Fi gure 5 and 
C7). Notably, even though two fewer constraints were 
imposed, they find a significantly smaller (2% versus 9%) 
constraint-satisfying volume. In large part, these differ- 
ences can be ascribed to incl uding systematic errors: in 
IO'Sha ughncss v~et al.l (|2005d ). preliminary data for the 
sparse and poor-quality NSNS(vw) sample was fit and 
applied without any account for fit-induced errors. How- 
ever, our calculation also differs at several fundamental 
levels from the original approach : (i) the space of mod- 
els studied is larger, covering more area in v\,V2 (see, 
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Fig. 3. — Probability distributions for merging (top) and wide 
(bottom) NS-NS formation rates, as predicted by population syn- 
thesis. The light shaded region shows the 95% confidence inter- 
val consistent with observations of binary pulsars, as described in 
Sec. [2] Because of sparse sampling and relatively poor data, our 
fit for the formation rate of visible, wide NS-NS binaries is signifi- 
cantly more uncertain than other fits. 



e.g., PSF2); (ii) the observed sample of merging NS-NS 
binaries is compared with the total merging NS-NS for- 
mation rate predicted from population synthesis, rather 
than with subset of merging NS-NS which undergo mass 
transfer; and critically (iii) the constraints are more nu- 
merous.^ 

As in lO'Shaughnessv et al.l (|2005ch . evaluating fits for 
other DCO merger rates on the constraint-satisfying 
model parameters provides revised estimates for various 
phenomena of interest, such as for the BH-NS and NS- 
NS merger rates. Figure [5] shows smoothed histograms 
of the BH-NS and NS-NS merger rates, both before and 
after observational constraints were applied. To be con- 
servative, these distributions also incorporate our best 
estimate for systematic errors associated with fitting to 
the data: rather than simply showing the histogram of 
merger rates that follows by evaluating the fits R(x) on 
many randomly poppulation synthesis parameters x, we 
smooth the resulting histogram by a gaussian with width 
J q (i.e., by the characteristic systematic error). For BH- 
NS and NS-NS binaries this characteristic error is sig- 
nificantly smaller than the characteristic width of the 
distribution. 

Finally, in Figure [7] we estimate the range of initial 
and advanced LIGO detection rates implied by these 
constrained results. This estimate combines our con- 
strained population synthesis results; our estimates for 
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Fig. 4. — Differential probability distributions dP^/dx defined so 
Pk(x) is the fraction of all models consistent with all DCO observa- 
tions and with x k < X. The top panel shows the distributions for 
the 3 kick-related parameters X3,X4,xs (dashed, solid, and dot- 
ted, respectively, corresponding to v\, V2, and s). The bottom 
panel shows the distributions for x\ (the mass-ratio distribution 
parameter r, appearing as the solid nearly constant curve), and 
the three binary-evolution-related parameters X2,xg,x^ (dashed, 
dotted, and solid, respectively, corresponding to to, «A, and f a ). 
The distribution of f a = £7 exhibits a strong peak somewhere 
between — O.f . Our choice of smoothing method causes all pro- 
jected distributions to appear to drop to zero at the boundaries, 
as it involves averaging with empty cells. 



systematic fitting error; the approximate LIGO range as 
a function of chirp mass M c : d = d Q (A / f c /1.2Af Q ) 5 ^ 6 
where d — 15^/3/2, 300 Mpc for initial and advanced 
LIGO, respectively; estimates for the average volume de- 
tected for a given species, as represented by ^M c 15 ^ 6 ^} = 

11.2,2.3Mp for BH- NS and NS-NS, respectively (cf. 
lO'Shaughnessv et all l2005d our average chirp masses 
are systematically higher because of our higher maxi- 
mum neutron star mass); and a homogeneous, Euclidean 
universe populated with Milky Way-equivalent galax- 
ies with den sity O.OlMpc" 3 , each forming stars at rate 
3.5M Q yr _1 (jKalogera et al.ll2001h . The high rate of dou- 
ble neutron star mergers [NS-NS (m)] demanded by ob- 
servations explains most of the differences between the 
relative likelihoods of various merger rates based on a 
priori (solid) and constrained (dashed) expectations. For 
example, the lowest aX < 0.1 eliminate too many NS-NS 
binaries to be consistent with the many observed Galac- 
tic NS-NS binaries. Additionally, the differences between 
this paper's results and those found in PSClare in part 
produced from the different assumptions (priors) about 



the relative likelihood of different population synthesis 
parameters: our distributions of constrained and uncon- 
strained merger rates differ from previous studies in large 
part because this study allows a broader range of super- 
novae kicks. 

We emphasize that the new formation and merger 
rates of double compact objects that we use are the re- 
sult of the revised popu lation synthesis code StarTrack 
(jBelczvnski et al.1 l2006bl ) . The underlying physics, and 
formation scenarios of BH- NS and NS-NS binaries w ill 
be discussed in full detail in iBelczvnski et all (|2007af ). 

6. CONCLUSIONS 

lO'Shaughnessv et al.l (|2005d) performed a proof-of- 
concept calculation to compare population synthesis sim- 
ulations with observations, interpret the resulting con- 
strained volume, and apply the constraint-satisfying pa- 
rameters to revise a priori predictions of astrophysi- 
cally critical rate results. In this paper, we demonstrate 
that even when systematic model fitting errors are ag- 
gressively overcompensated for, observational constraints 
still provide information about population synthesis pa- 
rameters. These results include both surprisingly good 
agreement with pulsar kick distributions and strong con- 
straints on at least one parameter involved in binary evo- 
lution. Since in almost all cases our systematic errors 
appear much smaller than those from observations, we 
are confident that now observational limitations, rather 
than computational ones, primarily limit our ability to 
constrain population synthesis results. In particular, 
we expect stringent tests of population synthesis models 
are possible, beginning by imposing more observational 
constraints than varied model parameters (seven in our 
case) . 

Our analysis remains predicated upon a correct iden- 
tification of all parameters and input physics critical 
for the formation of double compact objects. Ad- 
mittedly, our understanding of binary evolution con- 
tinues to evolve; however multiple population stud- 
ies from different groups over the years lead to very 
similar conclusions about the basic input physics that 
primarily aff ects the formation rates of double com- 
pact objects (IBelczvnski et al.ll2002t iHurlev et all 120021 : 
iPortegies Zwart fc Verbuntj 119961 iFrver et al.lll998ft . In 
this paper we use a relatively new version of the 
S tarTrack code th at is described in great detail 
irlBelczvnski et all (l2006bh . A forthcoming paper 
(|Belczvnski et al.ll2007af ) will discuss how these updates 
affect the evolutionary history of NS binaries in signifi- 
cantly greater detail. 

Our analysis also remains predicated upon a correct 
and complete interpretation of the observational sam- 
ple and associated systematic and model uncertainties. 
In some respects, however, our models for observations 
and their biases may be incomplete or not representa- 
tive. We effectively assume pulsar recycling is required 
to detect binary pulsars in our implicit hypothesis that 
only one pulsar in NS-NS systems, the one observed, was 
ever likely to have been detected. And finally we use a 
canonical value for pulsar beaming correction factor 
for several systems with known spins but unknown beam- 
ing geometry. 

And to be comprehensive regarding the factors ne- 
glected here, our population synthesis parameter study 
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has not varied over all plausible models. We have 
neither considered p urely polar supernova kicks (see 
iJohnston et alj I2005D nor changed the maximum NS 
mass from 2.5M nor even allowed for a heterogeneous 
birth population of single stars and binaries (a binary 
fraction of 100% is assumed in all models). And of course 
we continue to use a single fixed star formation rate for 
the Milky Way, rather than treat it as an independent 
uncertain but constrained parameter. 

However, all these limitations can be remedied by 
closer study. Thus, we expect that either (i) all DCO 
merger rates will become determined to within 0(50%) 
or better (based on the comparative accuracy to which 
we know NS-NS merger rates), or (ii) experimental data 



will conflict with the prevailing notion of binary popula- 
tion synthesis, revealing flaws and limitations in classical 
parametric models for binary stellar evolution. 
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Fig. 6. — A plot of the probability distribution for the BH-NS 
(top) and NS-NS (bottom) merger rates per Milky Way-equivalent 
galaxy, allowing for systematic errors in the BH-NS and NS-NS 
fits. The solid curves result from smoothing a histogram of re- 
sults from a random sample of population synthesis calculations 
with systematic errors included. The dashed curve shows the same 
results, assuming model parameters are restricted to those which 
satisfy all four DCO constraints (both WD-NS and both NS-NS 
constraints) described in Sections [2] and \3\ and summarized in Ta - 
ble |2] Compare with Figure 5 oflO'Shaughnc ssv et all l|2005d ). 
which uses substantially different priors. 



Fig. 7. — Range of expected detection rates for initial (top) and 
advanced (bottom) LIGO for BH-NS (solid) and NS-NS (dotted) 
binaries, based on models that satisfy obser vational constra i nts and 
under the assumptions outlined in O'Shaughncssy ct al. (2005c) 
(e.g., using a fixed chirp- mass spectrum and LIGO range for de- 
tection, as well as assuming Milky- Way like galaxies fill the universe 
with density O.OlMpc - 3 ). These estimates incorporate systematic 
errors due to the fitting procedure. 
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